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Abstract 



We have measured the coexistence curve of the binary liquid mixture n-heptane+nitrobenzene 
(C7Hi 6 + C6H5NO2) near its consolute point using an optical method. In particular, the crit- 
ical exponent (5 describing the coexistence curve was measured for this system. Previous 
experimental values of (5 for n-heptane+nitrobenzene were higher than the typical theoreti- 
cally calculated value, an unusual, although not unique, occurrence. In an effort to study this 
discrepancy, we have used an improved experimental apparatus for our measurements. We 
have taken special care to minimize temperature gradients and maximize the temperature 
stability of our thermal control system. We have also exploited features of a known optical 
method to analyze, thoroughly, sources of systematic errors. We measured an apparent value 
of (3 as 0.367 ± 0.006 and by a careful study of the known sources of error we find that they 
are not able to remove the discrepancy between the measured and the theoretical values 
of p. We also measured the critical temperature of the system at Tq = (291.80 ± 0.02) K 
(18.65 °C). 

Keywords: binary liquids; coexistence curve; critical exponents; critical phenomena; n- 
heptane+nitrobenzene; optical interferometry; mach-zehnder. 
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1 Introduction 



One of the most salient features of critical phenomena is the universal behaviour of the 
so-called critical exponents. In this work, we report results from an experiment, in which 
we measured the critical exponent (3 describing the coexistence curve of the binary liquid 
n-heptane+nitrobenzene. The former being lighter floats on the latter when the system is 
below its critical consolute temperature of about 18.9 °C pQ. Theoretical studies predict the 
values of the critical exponents to a very high degree of accuracy and several experimental 
approaches can verify those predictions (reviewed in [2j [3]). The critical phenomenon in 
question is the miscibility of the two liquids as the temperature is varied through the critical, 
or consolute, temperature To Typically, binary fluids exist as two immiscible phases below 
Tq, while above the critical temperature, the two liquids become completely miscible and 
only a single liquid phase exists [4]. If we call the two phases U and L (for 'upper' and 'lower', 
respectively) , the order parameter of choice for these systems is usually the difference between 
the concentration of one of the species (e.g., n-heptane) 0h in phase U and its concentration 
in the phase L [5]: A0 H = 0h,u — 0h,l- Introducing the reduced temperature t = (T c — T)/Tq, 
the shape of the coexistence curve very close to the critical point should be described by 
the simple scaling equation: A0 = B Q t", while further from the critical point correction-to- 
scaling terms are needed: A0 = B Q t^(l + Bit A + B 2 t 2A + •••), where the amplitudes B , 
£>i, B 2 , . . . are system-dependent coefficients to be determined, and A is the correction-to- 
scaling critical exponent [6]. The theoretical value of (3 is predicted to be approximately 
f3 = 0.326 ±0.002 [TJ. One way to test the validity of theoretical predictions on the exponent 
(3 is to measure the coexistence curve of a binary liquid mixture by optical interferometry, 
which is the method employed in this work [8]. 

The decision to measure the coexistence curve of a binary liquid mixture stems from the 
need to clarify discrepancies between the theoretical values of (3 and experimental results 
first obtained in the 1970s [91 [2]. It appears quite clearly from earlier literature that until 
the late 1970s or even early 1980s the typical accepted experimental values for the critical 
exponent (3 describing the order parameter dependence on temperature in binary mixtures 
was somewhat larger than the expected theoretical value for that exponent ([9J [21 EH EH])- 
In a compilation of experimental data on binary liquid mixtures, Stein and Allen [9] found 
that the data of the systems they analyzed were all consistent with the critical exponent 
(3 = 0.34 ± 0.01. The trend of other experimental studies of coexistence curves of binary 
liquid mixtures seemed instead to be confirming the renormalization group theoretical result 
that the critical exponent j3 should range between 0.325 < (3 < 0.327 [81 [T2"]. 

Following the relative uncertainty of the experimental results obtained on binary liquids, 
coupled with the fact that one of the systems displaying a greater than expected value of 
f3 was n-heptane+nitrobenzene, it seemed fit and interesting to perform a new study of the 
coexistence curve of that system, with a new and improved experimental apparatus. We 
believe that the thermal control system employed for this work allows greater temperature 
stability and uniformity. Moreover, the calibration of the thermometers was carried out more 
accurately and features of our optical technique were used to monitor more carefully the 
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Legend: 

1 He-Ne Laser 

2 Beam Expander and Spatial Filter 

3 Collimating Lens 

4 Polarizer 

5 Expanded Beam 
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Figure 1: Schematic diagram of the experimental optical setup used for the focal plane and 
image plane detection techniques. See figure legend for explanation of various parts of the 
setup. 

status of equilibrium (or lack thereof) of the system. All of these aspects of the experimental 
procedure were treated with greater care than in earlier studies of n-heptane+nitrobenzene, 
so as to obtain as accurate as possible a value of the critical exponent (3. 



2 Experimental features 

2.1 Optical observation of fluid critical phenomena 

The experimental method is based on the detection of variations that occur in the index of 
refraction of the liquid contained in a cell as its temperature varies. Our measurements are 
based on two optical detection techniques: the focal plane and the image plane techniques 
(Fig. [p. These techniques were already used for similar experiments in the past, and we 
refer to those studies for details [El [131 LT3] - In the focal plane technique, we exploit the 
refractive index gradient occurring as the sample's temperature is raised from below to 
above the critical temperature Tq. In the process, a fraction of the heptane diffuses into the 
nitrobenzene and vice versa. As a result, the refractive index profile, discontinuous below 
To, takes instead a sigmoidal shape (Fig. [2]). The refractive index profile of the sample is 
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Figure 2: Illustration of the behaviour of the refractive index of a binary mixture as function 
of cell height at T < T c (a), at T c (b), and T » T c . z marks the position of the meniscus 
between the two phases, below T c . 



probed using a coherent beam of light from a 632.8 nm He-Ne laser, with the beam expanded 
to approximately 25-mm in diameter and collimated. Light passing through the cell is bent 
due to the refractive index gradient in the sample. Rays encountering the sigmoidal profile 
at points with equal absolute value of the curvature are bent through equal angles. If the 
light is then focussed on a screen by a lens, a diffraction pattern appears in the focal plane of 
the lens. This diffraction pattern is recorded on photographic film. The number of fringes, 
N, observed in the diffraction pattern at each temperature is directly proportional to the 
difference in refractive index An between the two phases: An = n\j — n^. It can be shown 
that, to a good approximation under our experimental conditions [8j [15], An is proportional 
to the order parameter of the binary fluid system: An = kA<j) {k was previously measured 
as 6.045 [38]; we will discuss this issue further in section 14.41) . Therefore, a A0-versus- 
temperature graph represents the coexistence curve of the binary mixture and it is suitable 
for an indirect measurement of the exponent (3. 

In the image plane technique, the plane of the recording photographic film is placed on 
the image plane of the lens also used in the focal plane technique, and the whole refractive 
index profile of the sample is imaged on the film. As illustrated in Fig. [Tj this is achieved by 
placing the thermostat and sample in one of the beams of a Mach-Zehnder interferometer, 
while the other beam travels through air. As the refractive index varies considerably with 
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cell height, when the temperature of the sample is in the neighbourhood of the critical 
temperature, the bottom part of the sample has a larger optical thickness, hence the light of 
the sample beam that traverses the bottom part of the sample is slowed down more than the 
light passing through the top part. Therefore, when the waves at P and P' are recombined 
and imaged by means of the lens, horizontal interference fringes are observed at the image 
plane of the lens. The index of refraction as a function of height is mapped in this way. We 
used this technique to study the sample equilibration time at a set temperature. 

2.2 Sample preparation 

The liquids were purchased from Fisher Scientific with stated purities of 99% for nitrobenzene 
and 95% for n-heptane. They were distilled neat (the former under reduced pressure) and 
brought to an estimated purity of better than 99.9 %, as estimated by gas chromatography as 
well as by NMR. A mixture was prepared in the proportions of 49.6 wt% n-heptane and 50.4 
wt% nitrobenzene, as suggested by earlier studies [32] and the sample was prepared following 
a method devised and used by this laboratory in the past, with very good results [8]. 

Several different sample cells were prepared. The cells have a vertically elongated paral- 
lelepipedal shape, about 5 cm in height, 1 cm in width, and different depths (the dimension 
parallel to the He-Ne beam), depending on the light paths wanted for the measurements. In 
order to observe both the region far from the critical point and the near-critical region with 
relative ease, the sample cells were chosen to have light paths of 1, 2, 5, and 10 mm (the 
tolerance of the light paths as quoted by the manufacturer is 0.01 mm; the heights of the 
fluid in the 1-, 2-, 5-, and 10-mm samples were 34, 35, 37, and 38.5 mm, respectively). In 
our focal plane measurements, the number of fringes recorded at each temperature decreases 
as the critical temperature is approached, whereas the error made in counting the fringes is 
practically constant. Hence, in an effort to reduce the relative error on the fringe count for 
the data near the critical region, which are usually more significant, samples with longer light 
paths (5 and 10 mm) are used in that region. On the other hand, as data are taken farther 
and farther from the critical point, the number of fringes at each datum increases greatly 
for the same light path (~ 300 fringes at Tq — T ~ 0.1 K with the 10 mm cell), rendering 
the fringe counting operation more uncomfortable and prone to mistakes in counting. In the 
temperature region farther from Tq the 2- and 1-mm light path sample cells are employed. 

2.3 Thermal Control 

Thermal control of required stability and uniformity is achieved in two stages combined in 
a thermostat, whose detailed features are described in reference [17]. The external stage 
is regulated within ± 0.1 K, or better, at a temperature, T ext , such that Tq — T ext < 1 K, 
by circulating thermostated water in a water jacket system. The second, inner, stage is 
composed of a copper cylinder, to which heat is applied with heating foils, and of a cell 
holder, also a copper cylinder fitting snugly within the former. The inner cylinder and the 
cell holder are in turn placed inside the external stage. After the second heating stage, the 
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temperature can be controlled within about 5 x K. The temperature is measured by a 
quartz thermometer and by a set of thermistors embedded at various points in the copper 
cylinder. Both the quartz thermometer and the thermistors were calibrated using a triple 
point cell to reproduce the triple point of pure water. Tests were performed to ensure that 
adequate temperature uniformity and stability was achievable with this thermostat. By 
measurements done under different heating foil configurations, we estimated that thermal 
gradients along the height of the sample were less than 2 x K/m, which signifies a 
temperature difference of less than 6 x 10~ 5 K along the region of the meniscus, assuming it 
extends for a few mm in height. In terms of thermal stability, we monitored the temperature 
of the heater block, as measured by thermistors, for a time span of up to two weeks. In 
that time, the temperature of the experiment room increased by about 1.3°C, while the 
temperature measured by the thermistors remained stable to within less than 4 x 10~ 4 K. 



3 Results 

3.1 The critical exponent (3 

The collected coexistence curve data for the mixture n-heptane+nitrobenzene span 5 decades 
in terms of the reduced temperature t, from t ~ 4 x 10~ 7 to t ~ 4 x 10~ 2 . Data were 
taken so as to span the accessible range of temperatures of the mixture, which is around 
12.5°C and is bounded by the difference between the critical consolute temperature and the 
freezing temperature of n-heptane pQ. The data are reported in the online supplemental file 
accompanying this article. 

The coexistence curve data are shown in Fig. [3J, where they are expressed as (A^) 1 / - 326 
vs t. In this graph, we have plotted (A^) 1 / 0,326 as a function of the reduced temperature 
t, instead of the absolute temperature T, in order to be able to report on the same graph 
data from several experimental sequences with slightly different critical temperatures. We 
discuss below how the different critical temperatures arise. 

We anticipated in the introduction that in the vicinity of the critical temperature, the 
coexistence curve is supposed to be described by the simple scaling law: 

A0 = B iP (1) 

where t = 1 — T/Tq and (3 is the critical exponent. The "vicinity" of Tq, or t = — the so- 
called asymptotic region — has generally been found to be the region with t < 10~ 2 for binary 
liquid mixtures [2J. Translating this to the particular case of n-heptane+nitrobenzene, the 
region where the simple scaling law should be valid extends to about 3 degrees below Tq. It 
should be emphasized, however, that the asymptotic region is not known a priori, nor do we 
have a theoretical estimate of it. Moreover, in certain circumstances experiments on binary 
mixtures have suggested that the actual asymptotic region extends only to t ~ 10~ 3 from the 
critical temperature [12], [TBI [19] - In view of such findings, in the analysis of the present data 
a narrower asymptotic region < t < 10~ 3 was chosen to find the critical temperature and 
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Figure 3: The coexistence curve data for the binary mixture n-heptane+nitrobenzene ex- 
pressed as (A0//) 1 / - 326 as a function of the reduced temperature, where A(p H is the n-heptane 
volume fraction. The inset is an enlargement of the critical region. 



Fit 


Region 


B 


P 


x 2 


A 


t<i(r 3 


1.34 ±0.01 


(0.326) 


9.8 x 10~ 4 


B 


t<10~ 3 


1.91 ±0.04 


0.367 ±0.002 


1.3 x 10~ 4 



Table 1: Parameters values for a nonlinear least square best fit of A<f) = B^t 13 to the volume 
fraction data of Fig. [31 Quantities in brackets were held fixed during the fit. 



the order parameter critical exponent (3 from the measured data. Beyond the asymptotic 
region, correction-to-scaling terms become important: 

Acf> = Bot 13 (l ± B x t A + B 2 t 2A + ■■■) (2) 

where A is the correction-to-scaling critical exponent. 

Because it is physically inaccessible, the critical temperature Tq must inferred from the 
measured coexistence curve data. A first estimate of Tq is gathered by plotting the raw 
data as (A0) 1/,/3 versus T in the (expected) asymptotic temperature range, where such a 
plot is linear. A linear fit to the data intercepts the horizontal axis at Tq. This value of 
the critical temperature is then taken as an initial value for a nonlinear least square fit of 
([!)) to the coexistence curve data in the asymptotic region. In the fit, (3 and B are used 
as free parameters. (Tq is only allowed to vary within the reasonable range suggested by a 
careful examination of the coexistence curve near the critical point.) The best values of the 
parameters found by this nonlinear least square fit are given in Table [TJ As is apparent from 
the table, the fit obtained with (3 set at the theoretical value of (3 — 0.326 (fit A) does not 
seem to represent the data as well as the one where (3 is unconstrained. Instead, the value 
of the apparent exponent, which produces the best fit is (3 = 0.367 ± 0.006 (fit B). 
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Figure 4: Log-log plot of the order parameter A(p H versus the reduced temperature t. The 
slope of the data as t tends to zero corresponds to the apparent exponent (3. A line with 
slope 0.326 is also drawn for comparison. 
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Table 2: Parameters values for a nonlinear least square best fit of A0 = 
Bot 13 (l + 5 1 t A ±5 2 t 2A ), with t = (T - T c )/T, to the volume fraction data of Fig. [3j 
Quantities in brackets were held fixed at the shown value during the fit. 



Once an estimate of Tq has been obtained, a log-log plot of the data in the full temper- 
ature range is useful to see if any correction-to-scaling terms would be needed to interpret 
the data, and whether they should be positive or negative. Also, the slope of the data in 
the asymptotic region in the log- log plot corresponds to the critical exponent. A log- log 
plot for the data of Fig. [3] is reported in Fig. [H The slightly decreasing slope at higher 
values of t indicates that some correction to the simple scaling law are needed and that it 
will have to have a negative coefficient. This can be seen by performing a nonlinear least 
square fit of the scaling law (j2j) to the coexistence curve data of Fig. [3] in the whole range of 
temperatures studied, by using the best values found in the asymptotic range for the critical 
temperature Tq and for the critical exponent /3 and leaving the coefficient B\ free, while the 
correction-to-scaling exponent A is held fixed at its theoretical value of 0.54 [201 EI]- The 
line through the data in Fig. [3] corresponds to this fit and the parameters determined by the 
fit are reported in Table [2J fit C. 

It is evident from the best fit parameters of fits C and D that the scaling law with 
correction terms works better on the data when (3 is the value found by the simple scaling 
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Figure 5: Sensitive log- log plot, A(f> H versus t, of the coexistence data on n- 
heptane+nitrobenzene. The value used for the critical exponent is (3 = 0.326. A sensitivity 
scale is also drawn in the graph to indicate the slope the data would preferentially take in 
the asymptotic region, were they plotted with the value of the exponent corresponding to 
the indicated slope. 



law used in the asymptotic region (C) than when the theoretical value of (3 is imposed 
during the fit (D). When the exponent is again treated as a free parameter in the fit, with 
one correction term (fit E), the value of /3 agrees within error with the one found previously 
via the simple scaling law. 

A different way to perform the analysis of the data and arguably the ultimate verification 
of the importance of correction terms and of the goodness of the estimate of the measured f3 
is achieved by plotting the data in a manner that separates the contribution of the corrections 
terms from the rest of the terms in the relation (J2J). If no correction terms are needed to fit 
the data and the value of the critical exponent is "correct", a log-log plot of A^/t^ 3 vs t 
would distribute the data along a horizontal line (we call this a sensitive plot). Departures 
from a zero-slope line would then indicate either an "incorrect" value of f3 or the need of 
correction terms to fit the data, depending on the temperature range. Fig. [5] is a sensitive 
plot of the coexistence data collected during the experiment with n-heptane+nitrobenzene, 
where the theoretical value of (3 = 0.326 was used. It is evident once again that a value 
of (3 larger than the theoretical value is necessary to flatten the slope of the data in the 
asymptotic region. This type of graph can be used to perform a cross check on the values 
of the critical temperature Tq and (3 produced by a nonlinear least square fit to the raw 
coexistence data of Fig. [3l Starting from some good guesses for Tq and (3, one can then vary 
each of them individually step-by-step until the data of the sensitive plot lies on a horizontal 
line. The critical temperature will only affect the data very close to t — 0, while changes in f3 
will change the slope of the data on a wider range of t. The values of the critical temperature 
and the critical exponent produced by the experiment will then be the values that make, 
within experimental error, a zero-slope plot, in the asymptotic region. Moreover, if after 
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Table 3: Parameters values for a nonlinear least square best fit of A0 = 
Bot 13 (l + B x t A + B 2 t 2A ) (fit G), Acf) = Bot 13 (l + Bit A + B^t 1 - ) (fit H), and A0 = 

B t^ (l + Bit A + B'^t 213 ^ (fit I) to the volume fraction data of Fig. [31 Quantities in brackets 
were held fixed during the fit. 



this stepwise analysis the data appear distributed along different slopes at different ranges 
of t, this would be an indication that correction terms to the simple scaling law should be 
included in the fit. 

The modern theory of critical phenomena being as widely accepted and successful as 
it is, it seemed worthwhile to try some fits to the data by using the theoretical value of 
(3 and adding two correction terms in the scaling law (121). It has been suggested that dif- 
ferent exponents for the second correction term can be tried when analyzing coexistence 
curve data [231 121 [22]. Further fits to our coexistence data were then tried with the second 
correction terms being: B 2 t 2A , with A = 0.54, B' 2 t 1 ~ a , where a = 0.11 is the specific heat 
critical exponent, and B'^t 213 with (3 = 0.326. The rationale behind keeping the choice of the 
second term open is based on the experimental difficulty to distinguish between a correction 
exponent that is slightly larger (2 A) or smaller (1 — a or 2(3) than one. The results of these 
fits are reported in Tableland they appear to produce worse results than the experimentally 
measured apparent value of j3 — 0.367 ± 0.006. 

3.2 The critical temperature Tq 

Within each set of order parameter data the temperature was measured using the same 
thermistor, but different thermistors or other parts of the thermal control electronics were 
used in different sets. Moreover, as mentioned below in section 12. 2\ different samples were 
used to collect the body of data for this experiment. 

Systematic experimental studies [2U [25] on the effect of water and acetone impurities in 
different binary liquid samples (methanol±cyclohexane) show that a percentage volume of 
water of about 0.1 in the mixture would alter the critical temperature by almost ±4 K, while 
0.5% of acetone gives a increase of about 2 K. No regular variation of the critical temperature 
with the different samples was found and the different absolute values measured are all within 
a fraction of a percent of one another. While it is reasonable to think that both water and 
acetone impurities are present in the samples (the glass manifolds used to fill the sample cells 
are cleaned with both distilled water and acetone before making the samples), it is presumed 
that the amounts of those impurities do not differ much from one sample to the other. For 
these reasons, the actual critical consolute temperature of n-heptane±nitrobenzene varied 
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somewhat: while for each data set it can be determined to within less than 0.5 mK, its 
absolute value from our measurements can only be given as Tq = (291.80±0.02) K (18.65 °C). 

4 Discussion 

Due to the marked departure of our measured value of the critical exponent f3 from measure- 
ments on other systems as well as from the accepted theoretical value, we have carried out 
a careful analysis of sources of systematic errors, which may influence this type of measure- 
ments. We believe that the influence of other sources of error on the accuracy of our data 
was minimized by building an ad hoc improved thermostat and by being extremely careful 
in the conduction of the experiments and in the preparation of the samples. 

4.1 The effect of Earth's gravitational field 

Under the influence of Earth's gravitational field the refractive index vertical profile of a mix- 
ture below its consolute critical temperature is in principle distorted from a step-function 
shape like the one shown in Fig. [2^ (references [29j [26j [271 12H] address this issue in binary 
liquid mixtures in the one-phase region). It is conceivable that the gravity-modified profile 
would cause the appearance of diffraction fringes even when the system is at a temperature 
T < Tq, if the distortion from a pure step-like profile were appreciable. Although, no "un- 
wanted" fringes at temperatures below critical were ever detected in any of the experimental 
data sets, to study this point in more depth, the image plane method was employed because 
of its ability to map the index of refraction entire profile (see section |2~TI) . By counting the 
fringes as a function of their position on the film (the latter being related to the height in the 
sample cell), the profile can be plotted. Profile measurements were taken at a temperature 
above Tq, after the sample was shaken and therefore in a situation where a homogeneous 
index of refraction is expected, and below Tq roughly 50 hours after the temperature was 
lowered from above T c (Fig. Ej). To ensure that this method is sensitive enough to reproduce 
a sigmoidal profile if one is present, a measurement was made of the profile after the tem- 
perature was raised from below to above Tc, but without shaking the sample to homogenize 
the phase. In this situation a profile like Fig. [2b is expected, and is quite clearly measured 
by this method as is apparent from Fig. [7J From the measured profiles above and below the 
critical point, there seems to be no evidence of a dramatic deviation of the index of refraction 
from a homogeneous behaviour, both in the one-phase and the two-phase regions. 

Gravitational effect are supposed to be more evident the larger the density difference 
between the two species of the mixture [29]. The densities of n- heptane and nitrobenzene 
(Ph = 0.6837 and pn = 1.2037, at 20°C referred to the density of water at 4°C), while not 
as closely matched as those of other compounds [12], are however less mismatched than at 
least one other mixture that has revealed no influence of gravitational effect on analogous 
measurements [T8] . 
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Figure 6: Black diamonds: measurement of the refractive index profile at a temperature 
T > Tq after the sample was shaken to homogenize the phase. White diamonds: analogous 
measurement taken at a temperature T < Tc about 50 hours after the sample was cooled 
from approximately 8 x 10~ 4 K above to approximately 8 x 10 -4 below Tq. 
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Figure 7: Measurement of the refractive index profile at a temperature T > T c after the 
sample was heated from about 1 mK below to about 1 mK above critical, without shaking 
it. 
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4.2 Equilibration 

Closely tied to the effect of the terrestrial gravitational field, a very important issue to address 
is the equilibration time of a binary mixture. As a consequence of the divergences occurring 
at the critical point, the equilibration time also diverges as the set temperature T tends to 
Tq. Hence, the need to wait for an adequately long time after each temperature (<T C ) is set 
before heating the system through Tq to obtain a datum. (And consequently the need for 
a very stable, as well as accurate thermal control system.) Several other authors have put 
considerable effort in the study of this phenomenon [26j EH E3 EHJ ED] . From these studies, 
it appears that the most influential effect on equilibration is the sedimentation flux, which 
diverges as (T — T^) -0 ' 6 . Only once in the searched literature, however, did we find reference 
to the effect of this divergence on the two-phase region of a binary mixture, in a study by D. 
Beysens (within |30j). An order of magnitude calculation along the same lines suggests that in 
the asymptotic temperature range of our system the relaxation time due to sedimentation can 
reach several thousand hours (not unlike Beysens' findings), corroborating serious doubts on 
the feasibility of earring out meaningful critical phenomena studies of binary liquid mixtures 
on Earth. Clearly, this effect must have played a role in the experiments reporting results 
that are more compatible with the renormalization group theory predictions. 

We employed the image plane method (section 12. II) to estimate the equilibration time of 
the system, after a temperature change. This was done by monitoring the interference fringes 
at the image plane of the focussing lens with time as the film was moved in the camera. For 
our purposes, equilibrium was reached when the fringes appeared horizontal on the film. In 
the most significative measurement performed, the temperature was set at a value 7] slightly 
above Tq and held there for about two hours, it was then decreased to a temperature Tf 
below Tq, such that T\ — Tf <5 x 10~ 4 K. The system was kept at this temperature and the 
fringes recorded on film for several days. While the chart recorder track of the thermistors 
output showed that Tf was reached about 20 minutes after it was set, the fringes on the film 
do not appear to flatten to a horizontal slope until about 50 hours later. 

The coexistence curve data close to Tq were all taken after an equilibration time of about 
50 hours, while for the rest of the data the equilibrium time allowance was between 10 and 
20 hours. 

In principle, data taken before equilibrium is reached would cause one to measure a 
smaller number of fringes than the "true" number at a given temperature, hence a smaller 
value of the order parameter A0. This would, in turn, suggest that waiting a longer time 
would yield a more curved coexistence curve near Tc and hence a "more correct" value for 
P. 

4.3 Wetting 

There is some evidence that in binary mixtures one of the phases wets the walls of the 
cell containing the sample, as well as the other phase, sometimes surrounding the latter 
completely [31, 32J. The effect of this phenomenon in the experimental conditions of this 
work is that the effective thickness of the light path the laser traverses inside the cell is 
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different from the nominal thickness given by the manufacturer. Since the index of refraction, 
and hence the order parameter, is measured by the number of fringes detected at each 
temperature divided by the light path length, if the latter is made uncertain by the presence 
of a wetting film, the consequent refractive index measurements will be affected and rendered 
unreliable. 

By comparing refractive index measurements taken with cells of different thicknesses, the 
effect of a possible wetting layer can be monitored. Assuming the wetting layer that forms 
has the same thickness regardless of the nominal light path of the cell, if the data taken with 
different cells overlap within experimental error, we can neglect the influence of the wetting 
layer. The coexistence curve measurements were performed with different cell thicknesses to 
be able to map the whole range of t available for n-heptane+nitrobenzene with equal ease 
and accuracy. Those measurements as reported in Fig. H] and Fig. [5] indicate that the data 
indeed all seem to follow the same pattern. Within the accuracy of the measurements the 
data overlap, indicating that if a wetting film is present, it is, however, unmeasurable in 
these experiments. 

It is also worth emphasizing that there are indications from other observations [HI |2H [25] 
and theoretical predictions [33j that wetting behaviour occurs usually at temperatures several 
degrees below the critical temperature, a region of the coexistence curve with which one is 
less concerned when determining the critical exponent (3. 

4.4 Choice of order parameter 

As described above, we use index of refraction measurements to extract the volume frac- 
tion information to determine the coexistence curve of our mixture. In doing this three 
assumptions are made. The first is that the index of refraction does not present any singular 
behaviour at the critical point. There are theoretical studies [371 EH EHJ [36] and experimental 
observations dealing with this issue, both showing that any anomaly in the refractive index 
at the critical point is below 100 ppm, less than the resolution of these experiments. In 
the case of n-heptane+nitrobenzene, measurements were made in earlier studies to ascertain 
the goodness of the proportionality relationship between the mixture volume fraction and 
the order parameter [38]. The proportionality factor was found to be constant within about 
0.1% in the temperature range of the two-phase region of n-heptane+nitrobenzene. As such 
proportionality factor depends on the refractive index of the two phases of the mixture (see, 
for example, the reference within [15J), we feel safe in the assumption that in our system, 
too, any refractive index anomalies are likely not influential. 

The second assumption made is that of zero mixing volume. In other words, when the 
two individual species, H and N, are joined together to form the mixture it assumed that 
their volumes, Vn and Vn, simply add, while in general it is to be expected that V^xture = 
Vh + Vn ± |Ve|, where Ve is called the excess volume and it can be positive or negative. 
Experimental studies on two systems where Ve ^ have shown that this had no consequences 
on the measured critical exponent (3 [301 [35] . 

Thirdly, it is assumed that the Lorentz-Lorenz relation remains valid for binary mixtures 
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to the same extent as it valid for pure fluids, for which it holds within about 1%, at least in 
the neighbourhood of the critical point. According to an experiment aimed at determining 
the validity of this assumption [40l HI] H2], the Lorentz-Lorenz relation is verified within 
0.5%, when the volume loss upon mixing is considered. 

We have assumed that the above assumptions are valid for the results obtained with this 
system, although no measurements were carried out on it to verify such assumptions. 

5 Summary and conclusions 

Our measurement of the coexistence curve power law critical exponent for the mixture n- 
heptane+nitrobenzene as the temperature T — > T c yields an apparent value of (3 = 0.367 ± 
0.006, consistently higher than the accepted theoretical value of /theory — 0.326 ± 0.002 [7j. 
In an effort to find a reason for this discrepancy in some experimental flaw, the known 
potential sources of systematic error were carefully analyzed. In this process, the image 
plane optical technique was employed to measure the profile of the index of refraction of the 
binary liquid sample in a way that we had not tried before. The results are interesting in that 
the shape of the refractive index vertical gradient can be mapped directly by this technique. 
The relatively new use of this optical tool, as well as employment of an improved thermal 
control apparatus for the reported measurements have helped rule out that surface wetting 
by one of the phases, the influence of other effects such as impurities in the samples, and 
the exact definition of the order parameter could be responsible for the discrepancy between 
the measured and the theoretically predicted exponent (3. A quick literature survey revealed 
one other recent report of a similarly high value for (3 from an experiment [IB]- However, in 
that work, there appears no attempt to find a reasonable explanation for the result. 

Perhaps, two effects remain as candidate culprits for the high value of (3 we measured. 
On the one hand is the equilibration time issue. It is somewhat unclear (from this work 
and previous literature [26], [27J [28] [29] [BO]) whether and to what extent the influence of the 
gravitational field on the equilibration time of this mixture would affect the final outcome 
of our experiments. Estimates of equilibration times for a given system may vary by an 
order of magnitude (see, for example, remarks in [27] about observations obtained in |26j). 
Even in the least consevative case, those estimates would indicate that our equilibrium time 
allowances were too short and hence that some of our measurements, closer to the critical 
point, were performed before equilibrium was reached. Should that be indeed the case, while 
it would leave open a potential explanation for the discrepancy between our measured value of 
(3 and its widely accepted value, we would still be unable to explain why the fringes observed 
in the image plane measurements, aimed at establishing when equilibrium was reached, do 
become horizontal after a time supposedly shorter than the mixture's equilibration time 
(section 

On the other hand, it is conceivable that the asymptotic temperature range of this system 
be much closer to the critical point than expected for binary mixtures (or even for pure 
fluids). Our measurements in this case would suggest that the asymptotic region for this 
binary mixture would have to extend only to t < 10 -5 , a considerably narrower interval than 
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the more commonly observed t < 10 -3 . 
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